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1. Introduction 

A large number of materials are known which, over a certain temperature range, are well de- 
scribed by simple spin chain or ladder models [1-11]. In all these Mott insulators the superexchange 
constants are spatially very anisotropic so that a three-dimensional crystal effectively shows one- 
dimensional magnetic properties. In one of the best known spin-1/2 chain compounds Sr2Cu03, for 
example, the superexchange constant along the chain direction J ~ 2200 K whereas the magnetic 
couplings J± along the other directions are at least three orders of magnitude smaller [1,8]. In 
such a system one can therefore experimentally explore the physical properties of a spin chain over 
a very wide temperature range J_l ^ T < J. From a theoretical perspective this is very exciting 
because it allows to study experimentally many aspects of one-dimensional field theories [12-14]. 
In addition, the ideal Heisenberg spin-1/2 chain is integrable so that compounds well described 
by this model make it possible to experimentally address the question how physical properties, 
in particular, transport, are affected by a nearby integrable point [2-4]. Here the infinite set of 
constants of motion making the model integrable is expected to slow down the decay of current 
correlations - or even prevent them from decaying completely - leading to anomalous transport 
properties [15]. 

In real materials, however, we are always confronted with impurities and lattice imperfections 
which weaken or even completely destroy a superexchange bond between two spins. Because a weak- 
ening of a bond is a relevant perturbation in the renormalization group (RG) sense, we have to deal 
- at least at low temperatures - with finite chains with open boundary conditions (OBCs). Measure- 
ments on such systems then correspond to taking averages over ensembles of finite chain segments 
with lengths determined by a distribution function [16, 17]. Furthermore, the description by a one 
dimensional model will break down for temperatures T ^ J± where usually a three-dimensional 
magnetic order sets in. However, even for temperatures above this ordering temperature interchain 
couplings can have a significant effect which has to be taken into account. 

In the following, we will study the anisotropic spin-1/2 Heisenberg chain (XXZ model) with N 
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Figure 1. A non-magnetic impurity in (a) a single spin chain, (b) in a spin chain coupled to 
a neighboring chain by a ladder-like interchain coupling, and (c) in a spin chain coupled in a 
zigzag fashion to a neighboring chain. 
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Here J is the exchange constant and h the applied magnetic field. Although exchange anisotropics 
due to spin-orbit coupling are usually rather small so that experimentally only the isotropic case 
A = 1 is relevant, the additional parameter A is useful for the field theoretical calculations in 
section 3. There are two important consequences of the OBCs. First, expectation values of local 
operators become position dependent because translational invariance is broken. In the following 
we will in particular study the local susceptibility defined as 
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Second, there are well defined boundary contributions to all thermodynamic 
oo the total free energy F is, for example, given by 



F = iV/buik + Fb + 0{l/N) 



(1.3) 



where Fb is the boundary free energy. Similarly, one can define a boundary susceptibility xb [18- 
22]. These boundary or surface terms will be studied here as well. The local susceptibility defined 
in Eq. (1.2) can be related to the boundary susceptibility by 
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where Xbuik is the bulk susceptibility defined analogously to the bulk free energy in Eq. (1.3). 

For the interchain couplings we will consider two different cases relevant for many materials. 
One is a simple ladder-like antiferromagnetic coupling between neighboring chains as, for example, 
in Sr2Cu03 but also many other spin chain compounds. This case is shown in Fig. 1(b). The other 
is a zigzag ferromagnetic coupling between neighboring chains. This kind of interchain coupling is 
sketched in Fig. 1(c) and is relevant, for example, for SrCu02 [1]- 

To provide an intuitive picture we start with some numerical results for semi-infinite spin 
chains with geometries as shown in Fig. 1 obtained by the density matrix renormalization group 
applied to transfer matrices (TMRG) in section 2. In section 3 we will then present the field theory 
for finite Heisenberg chains with OBCs at finite temperatures. In section 4 we use the results 
obtained in the previous section to calculate the susceptibility as an ensemble average over finite 
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spin chain segments and show that the obtained results are in good agreement with experimental 
measurements. In section 5 we derive in a similar fashion the nuclear magnetic resonance (NMR) 
spectrum and show that it provides information about the interchain couplings. Finally, we give a 
brief summary and present some conclusions. 

2. Numerical results 

A method particularly suited to calculate the thermodynamic properties of one-dimensional 
systems is the density matrix renormalization group applied to transfer matrices (TMRG) [20- 
28]. To this end, the one-dimensional quantum system is mapped onto a two-dimensional classical 
system by a Trotter-Suzuki decomposition [29-31]. The additional dimension then corresponds to 
the inverse temperature /3. For the classical model a transfer matrix is defined which evolves along 
the spatial direction. Importantly, one can show that even for a critical system there is always 
a gap between the leading eigenvalue Aq and next-leading eigenvalues Aq of the transfer matrix 
T with = InjAo/Ao;! defining a correlation length. This makes it possible to perform the 
thermodynamic limit, i.e., system size N oo, exactly. With the TMRG one can treat impurity 
problems [21,22,32] as well as frustrated systems [33] making it the ideal numerical tool to study 
the problem considered here. 

We will start by investigating the local susceptibility as defined in Eq. (1.2) for a semi-infinite 
chain. By this we mean a chain which is infinitely long but has one end with OBCs. To obtain 
Xj we calculate the local magnetization (5?) for small magnetic fields h/J^ 10~^ and take a 
numerical derivative. Within the transfer matrix formalism the local magnetization is given by 

lim = mns^)Ti-'f\^v> . (2.1) 

Here T is a modified transfer matrix containing the broken bond, T{S^) is the transfer matrix 
with the operator S*^ included, and |\E'5j) ((^'lD the right (left) eigenstates belonging to the largest 
eigenvalue Aq, respectively. Far away from the boundary (Sj) becomes a constant, the bulk mag- 
netization 

m= hm lim (S^) = lim E„(^i|T(5-)T^-^|^".)(^2|f|n) 

Ao • ^^-^^ 

By taking again a numerical derivative with respect to a small magnetic field, we can obtain the 
bulk susceptibility Xbuik from (2.2). 

In Fig. 2 the susceptibility profile Xj ^ Xbuik for a single semi-infinite chain with A 1 is 
shown. As might be expected, the boundary induces Friedel-like oscillations which become larger 
with decreasing temperature. Interestingly, at low temperatures the oscillations first increase and 
reach a maximum, before decaying at large distances. This phenomenon has been first studied by 
Eggert and Affleck in [34] and we will rederive their field theory result as a special case of our more 
general considerations in section 3. 

Next, we consider the boundary susceptibility which we can easily obtain from the susceptibility 
profile using Eq. (1.4). The result is shown in Fig. 3. At low temperatures it was shown analytically 
that Xb ~ [Tln{To/T)]~^ with a known constant Tq [17-19]. The bulk susceptibility, on the other 
hand, behaves as Xbuik ^ const -I- ln~"'^(To/T) [12], i.e., it goes to finite value at T = with infinite 
slope. We will come back to this in section 3 but already notice that the numerical data are well 
described by the field theory. 

Finally, we also have calculated susceptibility profiles for the cases of coupled chains shown in 
Fig. 1(b) and (c). A ladder-Hke coupHng of neighboring chains is, for example, realized in Sr2Cu03 
[1]. In Fig. 4 it is shown that in this case an impurity in one chain also has a significant effect on a 
neighboring chain without impurities. The Friedel-like oscillations are refiected in the impurity-free 
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Figure 2. Xi " Xbuik for a semi-infinite chain with A = 1 at temperatures T/J = 
0.2, 0.1, 0.05, 0.025. The numerical TMRG data (closed symbols) are compared to the field 
theory formula (3.23) from section 5 (open symbols). 




Figure 3. The boundary susceptibility xs for A = 1 as obtained by TMRG (symbols) compared 
to the field theory result (3.20) (line) derived in section 3 which is valid at low temperatures. 
Inset: Numerical results for the bulk susceptibility Xbuik for A = 1 (symbols) compared to the 
field theory formula (3.19) (line). 



chain due to the interchain couplings. Clearly, the size of the reflected Friedel-like oscillations will 
depend on the ratio J±/T. However, this ratio is not the only relevant factor. The geometry of the 
interchain couplings plays an important role as well. In SrCu02 neighboring chains are coupled by a 
ferromagnetic zigzag-type coupHng. In this case an impurity in one chain leaves a neighboring chain 
almost unaffected even if the temperature T ~ Jj^ as is shown in Fig. 5. In section 5 we will show 
that these differences can be easily understood if one starts from the fleld theory results for a single 
chain derived in the next section and takes the interchain couplings into account perturbatively. 
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Figure 4. Xj — Xbu\k for A = 1, J_l = 0.03 J and T = 0.09 J and a ladder-like interchain coupling 
as shown in Fig. 1(b). One of the chains has a non-magnetic impurity at site j — (circles) 
whereas the other one (diamonds) is infinitely long and does not have any impurities. The 
numerical TMRG data (closed symbols) are compared to the field theory formula (5.4) from 
section 5 (open symbols). 




Figure 5. Xj ~ Xtuik for A = 1, J± — —0.1 J and T — 0.09 J with a zigzag interchain coupling 
as shown in Fig. 1(c). One of the chains has a non-magnetic impurity at site j = (circles) 
whereas the other one (diamonds) is infinitely long does not have any impurities. For this case 
field theory predicts that there are no oscillations, Xj ~ Xtuik ~ 0, in the infinitely long chain. 



3. Field theory for finite spin chains 



In the limit of low temperatures and large chain length, the XXZ model (1.1) can be represented 
by a field theory. The main step is a linearization of the dispersion around the two Fermi points. 
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A certain linear combination of particle-hole excitations around a Fermi point then constitutes a 
collective bosonic mode described by the Luttinger liquid Hamiltonian 



i-L+a nr pL+a 

H = - dx[ll^ + {d,^f]-hd— dxd,^ (3.1) 
^ Jo » ^^Jo 

where v is the spin velocity, L = Na, and a the lattice constant. The bosonic field $ obeys the 
standard commutation rule [$(a;), n(a;')] = iS{x — x') with 11 = v~^dt^. The Luttinger parameter 
K is & function of the anisotropy A and can be determined exactly by Bethe ansatz with K = 1 
at the isotropic point. The spin operators can be directly expressed in terms of the boson $, in 
particular, we have 




-l)^cosV2^$ (3.2) 

at zero magnetic field. Here c is an amplitude which can also be obtained exactly [35-37]. The 
separation of Sj into a uniform and a staggered part at low energies can be understood as follows. 
In the equivalent spinless fermion representation of the XXZ model, obtained by a Jordan- Wigner 
transformation, the S*^ operator becomes the density operator. Due to the linearization of the 
dispersion the electrons can live either at the same Fermi point and have therefore small momentum 
or they can be situated at different Fermi points in which case the associated momentum is 2kF- 
Zero magnetic field corresponds to half-filling for the spinless fermions thus 2kF = and we obtain 
the staggered contribution in (3.2). This also means that the local susceptibility defined in Eq. (1.2) 
can be separated at low temperatures into a uniform and a staggered part 

X. = X""' + {-lYxf ■ (3.3) 

In a bulk susceptibility measurement the staggered part does not contribute. This part is, however, 
measurable in probes of the local magnetism as for example in NMR. We will come back to this 
point in section 5. 

If one is only interested in correlation functions for infinite system size at finite temperatures 
or correlation functions for finite systems at zero temperature, one can make use of the conformal 
invariance of the field theory. It is then sufficient to calculate the correlation at zero temperature 
and infinite system size and use a conformal mapping from the complex plane onto a cylinder 
with the circumference corresponding either to inverse temperature or system size. Here, we are, 
however, interested in the thermodynamics of finite chains and we will therefore have to use an 
explicit mode expansion 

^(x=ja,t) = ^- + ^-5-,,^+E ^ ^ '' {^ ^-&.+e'-6t) (3.4) 

which incorporates the OBCs. Here 6„ is a bosonic annihilation operator. Eq. (3.4) is a discrete 
version of the mode expansions used in [13, 16] with x = ja becoming a continuous coordinate for 
a —^ 0, N CO with L = Na fixed. Using this mode expansion, the local observables respect the 
discrete lattice symmetry j ^ N +1 — j corresponding to a refiection at the central bond (site) for 
N even (odd), respectively. The sites and A'^ -I- 1 are added to model (1.1) and we demand that 
the spin density vanishes at these sites. Therefore the upper boundary for the integrals in (3.1) is 

L + a. The zero mode part (first two terms of Eq. (3.4)) fulfills « /o^^" = S^^^ 

and the oscillator part (last term of Eq. (3.4)) vanishes for j = Q, N + 1 as required. 

For the free boson model the uniform part of the susceptibility can easily be calculated and is 
given by [13] 



^2 

^° dh? 



1 Y.s.Sle-T^^^^"' 1 

JO = ~ 



M=0 



\ne [e'T^^f^ ,v}j (3.5) 
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Here 0{q,u) is the elliptic tlieta function of the third kind 6 — 9^{q,u) — J2^=-qo1" for 
integer Sz (even N) and of the second kind 9 = 92{q,u) = X^^.^.^ g^"+^/^''^e*(^"+^)" for half- 
integer Sz (odd A''). Note, that Xo"' has a simple scaling form as a function of NT. The lattice 
parameter a appears here due to the OBCs. It leads to a boundary correction which we will consider 
later. First, we set a = in (3.5) and concentrate the following limiting cases 

Tn [-wet] ^'^I'" ^ 0, at even 
xT = <! WN NT/v ^ 0, AT odd (3.6) 

The Curie- like divergence for N odd is caused by the degeneracy of the ground state S^^^ = ±1/2. 
At low temperatures the whole chain therefore behaves like a single spin. For TV even, on the other 
hand, the ground state is a singlet, S^^^ = 0. At low temperatures the chain becomes locked in this 
state leading to an exponentially small susceptibility. In the thermodynamic limit, NT/v — > oo, 
the susceptibility within the free boson approximation is just a constant. 

The staggered part of the susceptibility, xf ^ foi^ ^ finite chain with OBCs has been calculated 
in [37]. It is given by 

st^ c( TT \^ /%3K/2(g-f^) E,n^sin[27rmj7(iV+l)]e— 

Here t/(.t) is the Dedekind eta-function and 9i{u,q) the eUiptic theta-function of the first kind. 
The summation index m runs over all integers (half-integers) for N even (odd) , respectively. In the 
thermodynamic limit, N — > oo, we can simplify our result and obtain 

with X = ja. This agrees for the isotropic Heisenberg case, K = 1, with the result in [34]. The 
amplitude c, first introduced in Eq. (3.2), can be determined with the help of the Bethe ansatz 
along the lines of Ref. [36]. This leads to c = ^A^/2 with Az as given in Eq. (4.3) of [36]. The 
formulas (3.7) and (3.8) are therefore parameter free. 

To find the boundary contributions, as for example the boundary free energy defined in Eq. (1.3), 
one has to go beyond the free boson model (3.1). When deriving the low-energy effective theory 
for the XXZ model, one finds in addition to the free boson model (3.1) infinitely many irrelevant 
terms. These terms either stem from band curvature (corrections to the linear dispersions around 
the Fermi points) or from the interaction term. For A close to 1 the leading irrelevant term is given 
by 

rL+a 

5H = X dx cos{V8tt Kcj)). (3.9) 
Jo 

It is the bosonized version of Umklapp scattering where two left moving electrons get scattered 
to right movers or vice versa interchanging a reciprocal lattice vector. This term becomes relevant 
for A > 1 and is responsible for the opening of an excitation gap in this regime. For the isotropic 
case, A = 1, Umklapp scattering is marginally irrelevant and leads to important corrections to the 
results obtained for the free boson model. Due to the intcgrability of the XXZ model the amplitude 
A of the Umklapp term can be obtained exactly as well [38]. In Ref. [16, 17] the free energy and 
the susceptibility corrections to the free boson result to first order in the Umklapp scattering have 
been calculated. For the susceptibility the following correction was obtained 
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with 

9o{y,q) = ' .1 + — ^2 • (3-11) 

In addition, there is also the boundary correction related to the parameter a in (3.5). Expanding 
in this parameter to lowest order we find 

^Xr = ^52(e-^) (3.12) 

where 

g2(g)= ^^- 1% - \ ' i ■ (3.13) 

a plays the role of a lattice constant and its value can be determined for A < 1 [K > 1) by the 
Bethe ansatz and is given by [20, 21] 

a = 2-1/2 gjj^ [t^k/{4K - 4)] / cos [n/{4K - 4)] . (3.14) 

The uniform part of the susceptibility of a finite chain is therefore given by x""'(L, T) — Xo"'(a = 
0) + ^Xi^^ + '5x2"' • The corrections to the uniform zeroth order susceptibility Xo^^o, = 0) are 
important here because in the thermodynamic limit they give the boundary susceptibility Xb- In 
this limit 52 (e-'^''/(^^'^)) ^ K'^T'^L'^ /{2tt'^v'^) and therefore 

XB,i = lim Lxr = ^ ■ (3.15) 

This is just a constant contribution to the boundary susceptibility. Much more important is the 
boundary contribution stemming from i5x""'. Here we find 



Xs,2 = ^lirn^ixr' = -A(^-j B{K,1 - mi^^ ' '^i^'i^)] l^— j , (3.16) 

with B{x,y) = T{x)T{y)/T{x + y), ip'ix) — dip{x)/dx, and ip{^) being the digamma function. Note 
that for 1 < K < 3/2 (1/2 < A < 1) the boundary spin susceptibility Xb shows a divergent 
behavior ~ ijj'^-'^K temperature decreases. 

3.1. The isotropic point 

At the isotropic point Umklapp scattering becomes marginal and simple perturbation theory is 
no longer sufficient. In this case we have to replace the Umklapp scattering amplitude by a running 
coupling constant g{L,T) which obeys a known set of renormalization group equations [38] 

l/c, + ln(5)/2 = ln(vW7i"e^^''^^min[L,t;/T]) . (3.17) 

Here 7 « 0.577 is Euler's constant and for the isotropic case considered here, the spin velocity is 
V = J-k/2. The uniform susceptibility is then given by 

x""(iv,r) = xr + ^xr (3.18) 

wliorc K ^ \ + g{L,T)/2 in the exponentials of (3.5) and A g{L, T)/A in (3.10). In this case the 
parameter a in (3.5) is not determined by (3.14) and has to be used as a fitting parameter. 

In the thermodynamic limit we can again split the susceptibility (3.18) into a bulk and a 
boundary part. For the bulk susceptibility this yields the result first derived by Lukyanov [38] 

1 / g(T) 3ff3(T) 73 2\ , ^ 
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Here we have also added the correction from Umklapp scattering as well as a T^-term which 
stems from irrelevant operators with scaling dimension 4 describing band curvature. The running 
coupling constant g(T) is given by (3.17) with L = oo. In the inset of Fig. 3 this formula is displayed 
in comparison to the numerical results. For the boundary susceptibility, on the other hand, we find 

where we included second order corrections in g as derived in [17]. The comparison of this formula 
with a = 1.5 and numerical results is shown in the main panel of Fig. 3. 

Finally, also the result for the staggered part of the susceptibility given in Eq. (3.7) has to be 
modified at the isotropic point. Here we find [37] 



1 ( t: 773/2 (e-fr) ^^msin[27rmj7(A^+l)]e-'^"'"'(i-9)/(-^^) 



^'^ = ~(27r3^)V4T (.FTlj ^p^(^,e-35%) Y.m^-.v^^ii-a)/(,LT) 

(3.21) 

Now the running coupling constant g also depends on the distance of site j from the boundary and 
is given by 

l/g + ln(ff)/2 = In (min[Coa;, Co{L - x), v/^e^/^'+T/Tj) (3.22) 

where the constant Cq is not known and has to be used as a fitting parameter. Note, however, that 
for low temperatures and x, L ~ x 1 the value of Cq becomes irrelevant and our result for Xj 
therefore again parameter-free. If we consider the thermodynamic Hmit of Eq. (3.21) we find 

St ^ 23 y/* 1 X , , 

X?= -7- . . (3.23) 



A comparison of this formula with numerical data is shown in Fig. 2. The agreement is good but 
not perfect. The main problem here is that the temperature and the length scale set by the distance 
from the boundary are competing. The renormalization group equations, however, cannot be solved 
with both scales present. The formula (3.22) is derived in the limit when only one of those scales 
matter and the corrections can be significant if this is not the case. 



4. The averaged susceptibility and a comparison with experimental data 

In an actual crystal we have imperfections and impurities which limit the length of a spin chain 
segment. It is important to emphasize again that any weakening of a link is a relevant perturbation 
in the renormalization group sense. It is therefore expected to be a good approximation to assume 
that at low temperatures we have spin chain segments with open boundary conditions. As long as 
we do not know in detail how these defects are distributed it seems to be reasonable to assume 
a Poisson distribution, i.e., the probability of having an impurity at site j is not infiuenced by 
the other impurities. We might, however, expect this assumption to break down for large impurity 
concentrations where some sort of impurity order might set in. Using a Poisson distribution we have 
a normalized probability P{N) = p{l — p)^ to find a chain segment with N sites if the impurity 
concentration is p. Any measurement then corresponds to taking an average over this ensemble of 
chain segments. For the susceptibility, for example, we find 

Xp=p'E^(l-f')'^^W- (4-1) 

N 

Note that in a bulk measurement only the uniform part of the susceptibility contributes. The 
staggered part cancels out. Using the formula (3.18) we can immediately calculate the average 
susceptibility Xp- It is, however, very useful and instructive to derive a much simpler approximate 
formula. To this end, we notice that we have two different regimes for a finite spin chain. If 
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Figure 6. Measured susceptibility T(x — Xconst) for Sr2Cui-2,Pda:03+i with impurity concentra- 
tions X — 0%, 0.5%, 1%, 3% (crosses from bottom to top) from Ref. [5]. Here a constant Xconst) as 
given in table 1 has been subtracted from the experimental data. Subsequent curves are shifted 
by 5 X 10~^. For comparison theoretical results according to (4.2) with j — 1 are shown with 
an effective impurity concentration p as given in table 1. 



the temperature is larger than the finite size gap, T/J > l/N, we can approximate x(-^) ~ 
Xbuik + Xb/N. If, one the other hand, Tj J < 1/N we expect to see more or less the ground state 
properties of the finite chain. This means, according to Eq. (3.6), that there will be no contribution 
when the chain length is even while x(A^) ~ 1/{ATN) if the chain length is odd. We can therefore 
introduce a length Nc = ^J/T where this crossover occurs. Here 7 is a parameter which we expect 
to be of order 1. For the average susceptibility we can therefore write 



N odd 
p l-p 



N 



(4.2) 



4T2-P 



l-(l-p)^/^ +(l-p) 



l-p 



PI 

T 



Xbuik + PXB 



where Xbuik is given by Eq. (3.19) and xb by Eq. (3.20). 

An interesting experiment has been performed by Kojima et al. [5] where Palladium has been 
doped into the spin chain compound Sr2Cu03. The Pd ions then replace the Cu ions, act as non- 
magnetic impurities, and cut the chain into finite segments. A complication in the analysis of this 
experiment occurs because even the undoped Sr2Cu03 samples already have a substantial amount 
of chain breaks. It is believed that this is mainly a consequence of excess oxygen, i.e., we are in 
reality deahng with Sr2Cu03+5. In a simple picture, an excess oxygen ion pulls two electrons out 
of the copper chain. If these holes are relatively immobile this also corresponds to a chain break. In 
the comparison of the experimental data for Sr2C\ii-xP(^xO?,+s and formula (4.2) shown in Fig. 6 
we therefore use p as an effective impurity concentration incorporating both the chain breaks 
due to excess oxygen and due to the non-magnetic Pd ions. Furthermore, we have subtracted a 
constant contribution Xconst from the experimental data which is expected to be present due to 
core diamagnetism and Van Vleck paramagnetism. The values for p and Xconst which yield the 
best fit of the experimental data are shown in table 1. For Pd concentrations of a; = 0.5% and 
X — 1% the obtained values for p are consistent with the picture of having a certain amount of 
chain breaks due to excess oxygen already in the undoped compound. For x — 3%, however, this 
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Table 1. Concentration x of Pd ions in experiment compared to chain break concentration p and 
constant contribution Xconst yielding the best theoretical fit. The first line corresponds to the 
"as grown" sample of Sr2Cu03+5 from Ref. [1]. 



X (Exp.) 


p (Theory) 


Xconst [emu/mol] 


0.0 


0.006 


-7.42 X 10"^ 


0.005 


0.012 


-7.7 X 10"^ 


0.01 


0.014 


-7.5 X 10"^ 


0.03 


0.024 


-7.5 X 10"^ 



picture seems to fail. Reasons for this could be either on the experimental side (perhaps not all 
Pd ions really go into the sample, or some go in interstitially) or in the theoretical description. If, 
for example, the Pd ions tend to cluster above a certain concentration, then our assumption of a 
Poisson distribution becomes incorrect. 

5. The Knight shift and the role of interchain couplings 

As already mentioned, the staggered part of the susceptibility cannot be observed in a bulk 
measurement. This, however, is possible by NMR because here the resonance frequency gets shifted 
by the local magnetic field. This so-called Knight shift Kj can therefore be directly related to the 
local susceptibility. For a chain of length A'' it is given by 

<^ = (7e/7.)E^'"''xr (5-1) 

where 7e (7n) is the electron (nuclear) gyromagnctic ratio, respectively and Aj^y the liyperfinc 
coupling tensor. It is usually sufficient to take only AP and A^^ into account because of the short- 
range nature of the hyperfine interaction. 

To compare with experiment, we assume again a Poisson distribution of chain breaks so that 
the measured Knight shift is given by an average over all possible chain lengths. Furthermore, each 
site in a chain of length N gives a different Knight shift according to (5.1). If we assume that each 
of these Knight shifts has a Lorentzian lineshape with width T we find for the NMR spectrum 

Using the results for the uniform (3.18) and the staggered part (3.21) we can now immediately 
calculate the NMR spectrum for an ensemble of isotropic Heisenberg chain segments and compare 
to experimental data for Sr2Cu03 obtained by Takigawa et al. [6,7]. We use J = 2200 K as 
exchange constant in (1.1) and hyperfine coupling constants y4.°/(2fi,7„) —13 T, A°^/(2/i7„) fa 2 
T, and A^/(2/l7„) w 4 T [39]. Here the index denotes the magnetic field direction. We calculate the 
spectrum as a function of h = {\ + K)h^^^ and use h^^^ = 7.598 T. For the resonance field = 86 
MHz used in the experiment [6] this is consistent with /i°gg = v/^n where 7„ ~ 11.3 MHz/T [40]. 
A comparison of formula (5.2) for this set of material-dependent parameters with experiment is 
shown in Fig. 7. Here the impurity concentration p and the Lorentzian lincwidth T have been used 
as fitting parameters. Note that the sample used for the NMR experiment has been annealed thus 
dramatically reducing the amount of excess oxygen and the associated chain breaks compared to 
the sample used for the susceptibility measurement shown in Fig. 6. The theory predicts a central 
peak corresponding to the bulk susceptibility value and two shoulders with separation A/i (see curve 
'theory V in Fig. 7) which are caused by the maxima of the local susceptibility (see, for example. 
Fig. 2). Theoretically, we find that A/i ~ ^/vJT In^ ^ {v / T) , i.e., the separation of the shoulders 
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Figure 7. NMR spectrum for Sr2Cu03 at T = 30 K taken from [6]. In comparison the theoretical 
results for a single chain with a Poisson distribution of chain breaks (theory I) and for the case 
where also interchain couplings {J± = 5 K) are taken into account (theory II) are shown. In 
both cases p — 5 x 10"'', F = 4 x 10"*, and h°^^ = 7.598 T. The material-dependent parameters 
are given in the text. 



increases ^ 1/VT with decreasing temperature. This is in agreement with experimental findings 
[6]. The central peak, however, has a much more complex structure than predicted. The peak is 
split (feature 6h in Fig. 7) and has prominent shoulders (feature Dh in Fig. 7). These features 
are only observed in experiment at temperatures T < 60 K, whereas at higher temperatures the 
theoretical prediction agrees well with experimental data (not shown). 

To explain the additional features in the NMR spectra at T < 60 K we have to take the 
interchain coupling into account. Sr2Cu03 shows three-dimensional Neel ordering at T/v ^ 5 K. 
According to the usual chain mean field argument the interchain coupling should therefore also 
be of order J± 5 K. This estimate is consistent with band structure calculations [8]. As we 
have already demonstrated numerically in Fig. 4, a ladder-like interchain coupling as present in 
Sr2Cu03 leads to a reflection of susceptibility oscillations. Therefore another typical Knight shift 
is expected to be present in the NMR spectra, related to the maxima of the reflected oscillations. 
As long as J±/T ^ 1 we can calculate this effect perturbatively. Starting from the definition of 
the local susceptibility given in Eq. (1.2) we obtain to first order 

X^jA — {^j,l^k,l) {^tot,2^k,2) ■ (5-3) 

G"{j-k) Xk,2 

Here, the lower index 1 stands for the infinitely long chain without chain breaks whereas the lower 
index 2 denotes the chain with a chain break at site k — 0. Both, the two-point correlation Gl^{j—k) 
and the local susceptibility Xk,2 have a uniform and a staggered part in the low temperature limit. 

We therefore find ^ « x^f + {-^Vxfi^ with x^f = -J^xT'^'^^T"''''^ /T and 

4!^ = -jri-iy E xt2Gr'\j ~ k) . (5.4) 

k 

Here Eq. (3.8) has to be used for xt,2 while Gl^'''*{j-k) = {S^S^Y* = cV [-^ sinh(5f |j - fc|)] ^ is 
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the staggered part of the bulk two-point correlation function. In Fig. 4, formula (5.4) is compared 
with the numerical result and good agreement is found. 

If we take these reflections into account, then also the additional features in the NMR spectra 
are explained as shown in Fig. 7 (curve 'theory IF). In particular, the shoulders with separation 
Dh directly correspond to the maxima of the reflected susceptibility oscillations. The splitting 
of the peak Sh is of different origin. It is in fact not a splitting but rather a loss of intensity at 
the value which corresponds to the bulk susceptibility. The oscillations and reflected oscillations 
spread over the entire crystal at low temperatures so that there are simply no sites left which show 
bulk behavior. Rather interestingly, the peak usually associated with the bulk susceptibility value 
therefore turns into a dip at low temperatures due to the presence of chain breaks and interchain 
couplings. 

Finally, we also want to shed some light on the role played by the geometry of the interchain 
couplings. To this end, we consider an interchain coupling as shown in Fig. 1(c). Such a coupling 
is realized in SrCu02 with a ferromagnetic J± ^ [—0.1 J, —0.3 J] [1,8]. Using again flrst order 
perturbation theory we flnd in this case 

k " " 

G"(j-k) Xfc-i,2+Xfc,2 

Separating this into a uniform and a staggered part we flnd Xj i'*^^^ = ~2 Jj_x""''^'^''x2'"*^°V^ 

Xf^^ = -^{-ly E (Xm - xt-1,2) Gr''\j - k) . (5.6) 
'= ^ ^0 ' 

Therefore no reflections will occur in this case to flrst order in perturbation theory consistent with 
the numerical results shown in Fig. 5. This also means that an NMR spectrum for SrCu02 would 
not show any shoulders associated with the coupling to the nearest neighbor chain. Reflections in 
chains further away might, however, be still possible which then would again lead to additional 
structures in the NMR spectra at low temperatures. 

6. Summary and Conclusions 

We have investigated here how chain breaks and interchain couplings affect the physical prop- 
erties of spin chain compounds. A weak coupling of two Heisenberg chains is an irrelevant per- 
turbation in the renormalization group sense whereas the weakening of a bond in an otherwise 
homogenous chain is relevant. Open boundary conditions are therefore the stable flx point. Due 
to this reasoning it is expected that a wide class of perturbations like impurities or dislocations 
present in any real compound can be efi^ectively described at low energies as a chain break. An 
experimental measurement then corresponds to taking an average over an ensemble of flnite chain 
segments with open boundaries. 

By combining a low-energy effective fleld theory with Bethe ansatz we have derived parameter- 
free formulas for the thermodynamics of flnite spin-1/2 Heisenberg chains with open boundary 
conditions. Particular emphasis was put on a calculation of the susceptibility. Due to the broken 
translational invariance there exists a site-dependent staggered susceptibility in addition to the 
uniform site-independent part. However, even the uniform part is affected by the open boundary 
conditions in the sense that a surface contribution arises which is not present for periodic boundary 
conditions. 

We have shown that susceptibility measurements on Sr2Cu03 doped with non-magnetic Pd ions 
are well described by the theory presented here. One of the complications arising for this compound 
is, however, that even the undoped sample apparently already has a relatively large amount of chain 
breaks which are believed to be caused by excess oxygen. The impurity concentration used in the 
theory to flt the experimental data therefore differs signiflcantly from the nominal Pd concentration. 
The emerging picture nevertheless seems to be consistent - at least at low impurity concentrations 
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- with a fixed concentration of additional chain breaks aheady present in the undoped sample. It 
would certainly be of some value to obtain experimental data for this or some other spin-chain 
compound where the impurity concentration is well controlled and a direct comparison with theory 
therefore possible. 

Furthermore, we used the field theory to calculate NMR spectra. Importantly, the Knight 
shift is proportional to the local susceptibility so that the staggered site-dependent part of the 
susceptibility, which cannot be observed in a bulk susceptibility measurement, becomes observable. 
As has already been shown previously [6] , the staggered susceptibility leads to a broad background 
in the NMR spectra with edges caused by the maxima of the staggered susceptibility. For Sr2Cu03 
it has, however been found that the NMR spectra at low temperatures show puzzhng additional 
features which have been ascribed in [7] to a coupling to phonons. Here we have shown that these 
features quite naturally arise in a spin-only model if the known interchain couphngs are taken 
into account. It is also important to note that the geometry of the interchain couplings is crucial. 
As a specific example we considered in addition to the ladder-like interchain coupling relevant for 
Sr2Cu03 also the zigzag-like interchain coupling reahzed in SrCu02. In the latter case no additional 
structures in the NMR spectra related to this interchain couphng will occur. NMR experiments 
on spin chain compounds are therefore not only helpful to study impurity effects, they can also be 
used to investigate the geometry and strength of interchain couplings. 

Acknowledgments 

I want to thank all my collaborators on this and related topics, in particular, Ian Affleck, 
Michael Bortz, Sebastian Eggert, Andreas Kliimper, and Nicolas Laflorencie. 

References 

1. Motoyama N., Eisaki H., Uchida S., Phys. Rev. Lett., 1996, 76, 3212. 

2. Pratt F.L., Blundell S.J., Lancaster T., Baines C, Takagi S., Phys. Rev. Lett., 2006, 96, 247203. 

3. Thurber K.R., Hunt A.W., Imai T., Chou F.C., Phys. Rev. Lett., 2001, 87, 247202. 

4. Sologubenko A.V., Gianno K., Ott H.R., Vietkine A., Revcolevschi A., Phys. Rev. B, 2001, 64, 054412. 

5. Kojima K., et ai, Phys. Rev. B, 2004, 70, 094402. 

6. Takigawa M., Motoyama N., Eisaki H., Uchida S., Phys. Rev. B, 1997, 55, 14129. 

7. Boucher J. P., Takigawa M., Phys. Rev. B, 2000, 62, 367. 

8. Rosner H., Eschrig H., Hayn R., Drechsler S.L., Malek J., Phys. Rev. B, 1997, 56, 3402. 

9. Hase M., Phys. Rev. Lett., 1993, 70, 3651. 

10. Ma S., Broholm C, Reich D.H., Sternlieb B.J., Erwin R.W., Phys. Rev. Lett., 1992, 69, 3571. 

11. Riiegg C, et ai, Phys. Rev. Lett., 2008, 101, 247202. 

12. Eggert S., Affleck I., Takahashi M., Phys. Rev. Lett., 1994, 73, 332. 

13. Eggert S., Affleck L, Phys. Rev. B, 1992, 46, 10866. 

14. Giamarchi T., Quantum physics in One Dimension (Clarendon Press, Oxford, 2004). 

15. Zotos X., Naef F., Prelovsek P., Phys. Rev. B, 1997, 55, 11029. 

16. Sirker J., Laflorencie N., Fujimoto S., Eggert S., Affleck L, Phys. Rev. Lett., 2007, 98, 137205. 

17. Sirker J., Laflorencie N., Fujimoto S., Eggert S., Affleck I., J. Stat. Mech., 2008, P02015. 

18. Fujimoto S., Eggert S., Phys. Rev. Lett., 2004, 92, 037206. 

19. Furusaki A., Hikihara T., Phys. Rev. B, 2004, 69, 094429. 

20. Bortz M., Sirker J., J. Phys. A: Math. Gen., 2005, 38, 5957. 

21. Sirker ,J., Bortz M., J. Stat. Mech., 2006, P01007. 

22. Sirker J., Bortz M., Phys. Rev. B, 2006, 73, 1. 

23. Bursill R., et ai, J. Phys: Cond. Mat., 1995, 7, 8605. 

24. Wang X., Xiang T., Phys. Rev. B, 1997, 56, 5061. 

25. Shibata N., J. Phys. Soc. Jpn., 1997, 66, 2221. 

26. Sirker J., Kliimper A., Phys. Rev. B, 2002, 66, 245102. 

27. Sirker J., Kliimper A., Europhys. Lett., 2002, 60, 262. 

28. Sirker J., Kliimper A., Phys. Rev. B, 2005, 71, 241101(R). 

29. Trotter H.F., Proc. Amer. Math. Soc, 1959, 10, 545. 

30. Suzuki M., Commun. Math. Phys., 1976, 51, 183. 



14 



Spin-1/2 chain compounds 



31. Suzuki M., Phys. Rev. B, 1985, 31, 2957. 

32. Eggert S., Rommer S., Phys. Rev. Lett., 1998, 81, 1690. 

33. Kliimper A., Raupach R., Schonfeld F., Phys. Rev. B, 1999, 59, 3612. 

34. Eggert S., Affleck L, Phys. Rev. Lett., 1995, 75, 934. 

35. Affleck L, ,J. Phys. A, 1998, 31, 4573. 

36. Lukyanov S., Terras V., NucL Phys. B, 2003, 654, 323. 

37. Sirker J., Lafiorencie N., arXiv:081 1.2106, 2008. 

38. Lukyanov S., NucL Phys. B, 1998, 522, 533. 

39. Monien H., Pines D., Takigawa M., Phys. Rev. B, 1990, 43, 258. 

40. Abragam A., Bleaney B., Electron Paramagnetic Resonance of Transition Ions (Clarendon Press, 
Oxford 1970). 



15 



